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Abstract. We study the phase— space of FLRW models derived from Scalar- 
Tensor Gravity where the non-minimal coupling is F(<f>) = £<f> 2 and the effective 
potential is V(<p) = \<f> n . Our analysis allows to unfold many feature of the 
cosmology of this class of theories. For example, the evolution mechanism towards 
states indistinguishable from GR is recovered and proved to depend critically on 
the form of the potential V(<f>). Also, transient almost-Friedmann phases evolving 
towards accelerated expansion and unstable inflationary phases evolving towards 
stable ones are found. Some of our results are shown to hold also for the String- 
Dilaton action. 



PACS numbers: 98.80.JK, 04.50.+h, 05.45.-a 

1. Introduction 

Scalar-Tensor Gravity (STG) has become a major area of activity in the last thirty 
years. This is due mainly to the discovery of several shortcomings of General 
Relativity (GR) in cosmology and quantum field theory. In fact, the presence of a 
big bang singularity, the flatness and horizon problems, led to the realization that 
the "standard" model for cosmology [1], based on GR and the standard model of 
particle physics, is not able to describe the universe in extreme regimes. On the other 
hand, GR is a classical theory which has so far eluded every attempt to be quantized, 
making the description of gravity on microscopic scales very difficult. This has led 
many authors to believe that the standard Einstein scheme is inadequate to fully 
describe the gravitational interaction at a fundamental level [2] and has stimulated 
the proposal and the study of alternative theories of gravity that might be able to 
solve these problems and at the same time reduce to GR in the weak field regime. 

One of the most fruitful theories of this type is STG, in which scalar field(s) 
that are non-minimally coupled to gravity are introduced. These type of theories 
were proposed nearly half a century ago by Jordan [3] and later refined by Brans and 
Dicke [4] . The original motivation behind Brans-Dicke theory (BD) came from Mach's 
principle, but over the years, BD and STG have gained interest in a wide number of 
scenarios. For example, in unification schemes such as superstrings, supergravity or 
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grand unified theories, the one-loop or higher-loop corrections in the high-curvature 
regime take the form of non-minimal couplings to the geometry or higher-order 
curvature invariants in their low energy effective actions [6]. 

In cosmology, STG acquired considerable interest because they introduce 
naturally scalar fields and scalar fields which are capable, among other things, of 
giving rise to inflationary behaviour [7]. In GR the introduction of this type of 
field has the drawback of raising the issue of explaining its origin. Instead, in STG 
this problem finds a natural solution, because the scalar field can be considered an 
additional degree of freedom of the gravitational interaction. For these reasons, among 
others, inflationary models based on STG have been widely studied [8,9]. 

Recently, STG has also been used to model dark energy, because scalar fields are 
also the natural candidates for phantom and quintessence fields (see [10] for a review). 
This suggests that both inflation and dark energy could be the result of the action of 
the same scalar field. STG offers, in this sense, the ideal framework to implement this 
intriguing idea and might even allow us to overcome the problems due to the energy 
scale difference between the dark energy scalar fields and the inflaton [11]. 

In this paper, we consider the problem of determining the global dynamics of 
Fricdmann-Lcmaitre-Robertson- Walker (FLRW) cosmologies of STG. We will study 
a generic class of STG theories, where quadratic non-minimal couplings to gravity 
and self-interaction power-law potentials are assumed. This class of models is strictly 
related to the String Dilaton action and naturally exhibit duality in the cosmological 
solutions [12]. Furthermore, they can be obtained from generic non-minimally coupled 
scalar-tensor Lagrangians if Noether symmetries are found in the dynamics [13,14]. 
Several exact solutions of these models have been found, but the stability and global 
behaviour is still not well understood. Our aim is to give a full description of the global 
dynamics of this class of STG and determine if cosmic histories are possible that (i) 
present a transient matter-dominated Fricdmann phase and then evolve towards an 
accelerated (dark energy-dominated or ACDM) regime or (ii) present a first unstable 
inflationary phase and a second inflationary attractor. 

The dynamical systems approach has been already used with great success in 
BD [16-20], in specific models of STG [21, 26-30] and in other contexts [31-33] 
allowing the derivation of new exact solutions and making possible a qualitative 
global description of the dynamics. The geometric structure and dimensionality of 
the phase-space of general couplings and potentials has recently been investigated by 
Faraoni [34]. 

The outline of this paper is as follows. In section 2, we derive the field equations 
and discuss the non minimal coupling and the self— interacting potential. The 
cosmological equations for a FLRW metric are derived in section 3. We convert them 
into an autonomous dynamical system in order to pursue a phase-space investigation. 
Sections 4 and 5 are respectively devoted to the analysis of this system in the absence 
and presence of a perfect fluid source. We study the stability of the fixed points with 
the aim of identifying physically interesting cosmological behaviours. In section 6, the 
results are discussed and some conclusions are drawn. 

The following conventions are adopted: the metric signature is (H ); Latin 

indices run from to 3; V is the covariant derivative; physical units c = 8irG = 1 are 
used. 
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2. The field equations 

In four dimensions, a general action in which gravity is non-minimally coupled to a 
scalar field 4>, reads [14]: 

A = [ dx^J—g [F(cf>)R + \g ab V a ^ V fc </> - V{4>) + C M ] , (1) 

where F(cj>) is a generic coupling, V(<p) is the self-interaction potential and Cm is the 
matter Lagrangian. 

By varying the action ([I]) with respect to the metric g a b, we obtain the field 
equations 

F{4>)R ab - \F{4>)Rg ab + ±V Q V 6 - \g ab V c <f> V c 

+ \g ab V{<P) - V a V b F(cp) + g ab UF{4>) = -T% (2) 

and the variation with respect to <f> gives the Klein-Gordon equation 

□0 - RF'(4>) + V'(4>) = 0, (3) 

where primes denote differentiation with respect to <fi. The field equations ^ can be 
recast in the standard form 

Gab = Rah ~ \R9ab = T a l , 

where the effective stress-energy momentum tensor is given by 

Kb = [-|V O V ^ + i9a fc V c V c - \g ab V{$) 

+ V a V b F(0) - g ab UF{4>) - T%] . (4) 

This form of the equations shows clearly one of the most interesting features of Scalar 
Tensor Gravity: the Newtonian gravitational constant Gjv, intended as the coupling 
of gravity with standard matter, has to be replaced by an effective gravitational 
"constant" G e // that depends on the non-minimal coupling F (and, as a consequence, 
on <t>) and varies in time. 

The Bianchi identities G a b - a = give the conservation laws for both the matter 
and the scalar field. As a general result [14], it is possible to show that the conservation 
law for the scalar field is the Klein-Gordon equation. 

The action |T]) is very general and contains several interesting physical cases. For 
example, considering the transformation 

4> cxp[-V], F{4>) - i ex p[-2# V{4>) -> ^Aexp[-2# (5) 

which specifies a particular form of coupling and potential, leads to the 4D-String- 
Dilaton Effective Action 

A = J dx i ^g~e~^ [R + g ab V a i' V b ip - A] , (6) 

where A is the string charge. In this context, such a theory is nothing else but a 
particular STG [35]. This means that the considerations and results for the action 
([1]), and the related dynamics, also hold for string-dilaton cosmology. On the other 
hand, the set of transformations 
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give rise to the action 



A= dx 



(8) 



which is nothing else but the BD actioilj] [4]. In addition to the ones above other 
interesting kind of transformations are possible [5] . 

In this paper, we will focus on models where the non-minimal coupling has 
the form F(<f) = £</> 2 and the self-interaction potential is taken to be an arbitrary 
power law of the form V((f>) = X<p n . This choice is general and motivated by several 
mathematical and physical reasons. In particular, beside the string-dilaton and BD 
actions, several effective quantum field theories, in low energy physics, can be related 
to such couplings and self interacting potentials [2,36]. Furthermore, this coupling 
and potential satisfy the requirement of Noether symmetries for the Lagrangian in 
([lj, giving rise to general exact solutions of physical interest [13]. 

A final remark concerns the parameters of the theory Very different models can be 
parameterized by the set (£, A, n), but not all the combinations of these parameters are 
necessarily physical. For example, attractive gravity is achieved for £ < and A > 0, 
although physically interesting situation can also be achieved for £ > 0, A < [37]. In 
this paper, we will consider only the physical case £ < 0, A > 0. The discussion of the 
other interesting combinations of parameters can be found in [38] . 



3. The FLRW dynamical system 

In order to analyze the phase-space of the Scalar Tensor FLRW cosmologies, the field 
equations need to be recast in a dynamical system form. In the FLRW metric and with 
our choice of coupling and potential the the above Einstein-Klein-Gordon equations 
reduce to: 



n k <b A „ o 1 <b a , . 

tp <p \a a z a z J 

where the first is the cosmological equation for the acceleration, the second is the 
Hamiltonian constraint, i.e. the {0, 0} equation, and the third is the Klein-Gordon 
equation. H = a/a is the Hubble parameter and k is the spatial curvature constant. 
We have also assumed standard matter to be a perfect fluid with a barotropic index 
w, so that the conservation equation yields 

(t = -3Hn(l + w), (12) 

where \i is the the matter-energy density. In what follows, since we want to stress the 
role of a non minimally coupled scalar field in the modelling of dark energy, we will 
consider only < w < 1. Other values of these parameters, which can be associated 
with more exotic forms of matter energy densities, although interesting, will not be 
considered here. 



If. To be precise, the proper BD action is exactly recovered only for w = constant. 
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The equations above can be converted into an autonomous system of first-order 
differential equations by defining the following set of expansion normalised variables: 

j,n-2 



<t>H'- 



y = 

K 



6£tf 2 ' 
k 



l H 2- 



from which we obtain 
f 



12f- f 

2y 



12£ (1 + K + Kx) + 2(1 + 6£)x + (5 - 24£)s 



1-AU 3 



66 

z 



+6£(n - 2)y - (1 - 6n£)xy + 6£(1 + 3w)z + -(1 + 3w)a;z 
[24£ - f + 12£K + Ax + (12£ - l)(n - 2)a: 



1 



1- A .T 2 -(f-67^)y + 7r( 1 + 3w ) 



12£ 
-2(l- 



- [7(1 - 120 - 2 - 3w(l - 120 + + 6(1 - 4C)a 



(l - ±"jx 2 - 2(1 - 6<)y + (1 + Zw)z 



2K 



12£- 1 



12£(1 + K) + 2x - ( f - ^ ) x 2 - (1 - 6n£)» + 7; (1 + 3w> 



(13) 



(14) 



where primes denote derivatives with respect to a new evolution variable r = In a and 
the dynamical variables arc constrained by 

1 + 2x + -jt^x 2 + y + K - z = 0. (15) 

The associated phase-space is 4-dimcnsional and the evolution is constrained by (If 5|) . 
The task is now to study the structure of such a space: this means finding the stability 
of the fixed points, and then to analyze the evolution of trajectories [15]. We will 
consider two cases: the vacuum case (p — 0) and the matter case (/1 7^ 0). 



4. The vacuum case 



When we consider the vacuum case (// = 0), the set of dynamical equations (|T4l 
reduces to 

m - 1 

+6£(n- 2)y- (1 - 6n£)xy] 



12^ (1 + K + Kx) + 2(1 + 60a; + (5 - 240x 2 - ( 1 - ^ ) x 



y 



2/y 



12£-1 



[24£ - 1 + 12£K + 4x+ (12C - l)(n - 2)ai 



- l-±)x 



(6<- l)y 



2A" 



12£(1 + A') + 2x - (l - x 2 + (6< - l)y 



12£-1 

with the constraint equation given by 

1 + 2x + -^x 2 + y + K = Q. 
In this case, the phase-space is 3-dimcnsional 



(16) 



(17) 
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4-1- Finite analysis 

We can further simplify the system (fT!)|) by implementing the constraint (JT7J): 
x' = -2x - 4x 2 - ±x 3 1 



+ I ^ rT [6^(n-4) + (6e(n-2)-l) a; ], (18) 



y = y 



V > 3i U £ _ 1 

From these equations, it is clear that the x axis (y = 0), characterized by the absence 
of a potential for the scalar field, is an invariant submanifold. This tells us that there 
is no orbit for which the potential of the scalar field can become exactly zero and that 
if the potential is initially set to zero, it will remain zero. 

The fixed points can be obtained by setting x' = and y' = 0. For the system (|18[l 
we obtain five fixed points (see TablcQ]). The coordinate of the point A is independent 
of £ and n and the ones of B and C are independent of n. The point T> is a finite fixed 
point for n^2 and £ a finite fixed point for £ ^ l/2(n + 2). 



Merging occurs for the points £ and B. and £ and C for n = 4 =F 3 \A( 12 £ 1 ) 



respectively. The point £ also merges with V for n = 4 ± y 12 | 1 . 

All the fixed points except A and 2? are associated with flat spatial geometry. 
Point A is associated with an open spatial geometry and for T> the sign of the space 

curvature depends on £ and n: K-p is positive for n ^ 2, 4— y^ 12 ^ -1 < n < 4+^/^2^! 

and negative otherwise. 

The stability of the fixed points can be determined by evaluating the eigenvalues 
of the Jacobian matrix associated with the system (fT5)) (sec Table [3]) , as prescribed 
by the Hartman-Grobman Theorem [39] |§[ 

The fixed point A is a saddle for every value of the parameters £ and n. The 
point B can either be a stable node or a saddle node whereas C is either an unstable 
node or a saddle node, depending on the values of £ and n. The eigenvalues of T> and 
£ are both dependent on £ and n so that the stability varies over the different ranges 
of £ and n. The stability of these fixed points have been summarized in Table [4j 

The coordinates of the fixed points may be used to determine exact cosmological 
solutions at the fixed points themselves. In fact, when evaluated at these points, the 
Friedmann and the Klein-Gordon equations can be written as: 



2 



H 

H = , a 



1 2x ? 



| jS H 2 , 13 = 3x 4 + x\ + 6 -B^yu (20) 

where (x.^yi) represent the coordinates of the fixed points. Integrating (|19p and 
substituting into (|2U)) gives 

f + (d^ = ' ^ < 21 » 

which has a Cauchy-Eulcr form. If the terms a and f3 arc different from zero, equations 
(|19p and (f2"Tj) can be easily integrated, giving 

a = a (t-t ) a , (22) 

§ The values of the parameter for which the eigenvalues are zero are bifurcations for the dynamical 
system. In this paper we will not give an analysis of the bifurcations referring the reader to the 
specific literature for more details. 
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(t - t ) 1/2 [0o (t - t ) m + 0i (t - t )- m \ ,, if a 2 /? < \ , 

(t - t ) 1/2 (0o + 0i lnt) , if a 2 /3 = | , ( 23 ) 

(i - io) 1 ^ 2 (0o sin [mln (i - to)] + 0i cos [mln (f — i )]) , if a 2 /? > J , 



where 



m = -Aa 2 (5. (24) 

In the case of point A, we have a Milne evolution and a constant scalar field 

a A = 1, 0.4 = O . (25) 

The above solution (as well as some others that will follow) is particularly interesting 
because the scalar field is constant. This implies that A represents a state in which 
G e ff is constant and the potential of the scalar field acts as a cosmological constant. 
In other words, at this point, scalar tensor gravity is indistinguishable from standard 
GR plus cosmological constant. 

For the point B the solutions are given by (|22[) and (f2"3")l with 

1 3 
ajs = = , me = = , (26) 

3(l-80-V3f(12f-l) 6(1-80-8^(12^-1) 
and for point C 

a c = == j mc — ; • (27) 

3(l-8£) + V3£(12£-l) 6(1-80 + 8^(12^-1) 

Note that, for ^ < 0, the value of a for the above two solutions is always positive and 

less than 1 i.e. these two solution always represent two Fricdmann-likc solutions. In 

addition, B represents a solution in which the scalar field is growing, while at C, is 

dissipating. If n ^ 2, the point V is associated with 

n + 2 

a v = 1, m v = ^- Yy (28) 

which represents another Milne solution, while the scalar field is decreasing for n > 2 
and increasing for n > 2. It is interesting that, unlike the Milne solution in GR, this 
linear solution for a is not necessarily a spatially hyperbolic one. The constant ao can 
be related to the parameters £ and n for non-flat solutions: 

2 _ k ( n - 2 ) 2 £ 

fl ° " -l + [4 + n (n-8)]f (29) 

When n = 2, T> becomes an asymptotic fixed point and merges with X>oo (its solution 
will be presented within the asymptotic analysis). Finally, for point £ we have 

2(n + 2)g-l » + 2 

(n-4)(n-2)£ 2(n-2) 

for all n 7^ 2, 4. This solution represent an expansion for < n < 2 and n > 4 and 
was already found in other contexts [21-25]. In particular, in [25] it is shown that for 
2 < n < 4 the scale factor associated to this point evolves towards a superinflating state 
(also called "Big Rip" singularity) without including any exotic feature like ghosts or 
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Table 1. The coordinates and scale factor solutions of the fixed points for the 
vacuum case. We only show the exponent a of the solutions (1221 . 



Point Coordinates (x, y) 



K 



A (0, 0) 

B (2(-6£-V3£(12£-l)), 

C (2(-6£+V3£(12£-l)), 

p ^ 2_ 2-24g 



/ 2(n 
U"2( 



-2' 3?(n-2) 2 

4)5 (12g-l)[3+ 



-2(n+2)£ > 



3(l-2(n+2)5) a 



-l-[4+nQt-8)]£ 



3(l-8«)-V3$(12e-l) 

1 

3(l-8C)+4^/3?(12£-l) 



fc(Tt-2) 2 g 



l+[4+n(n-8)]f y 

2(»+2)g-i 4 

(n-4)(n-2)£ ' " r * 

a = ao, n = 2 

a = aoe C(t_to) , n = 4 



non standard fluids. When n = 2 and n = 4. using the cosmological equations wc 
obtain the solutions 

a = a , cf> = (f) Q (A = 0) , (31) 

a = a e c (*-*o) ) = ^ o = ±v /lMiI 5 (32) 

which represent a static universe and a de Sitter evolution respectively. In both these 
cases an effective cosmological constant is present, whose value depends on the effective 
gravitational constant (via <f>o) and the coupling constant A of the self- interaction of 
the scalar field. Again, since the scalar field is constant in these cases, these solution 
are indistinguishable from the GR solutions. The difference with point A is that the 
solution does not occur in "pure" GR and can also be stable. This is particularly 
interesting in the second case (n = 4) in which a dc-Sittcr solution able to mimic 
an inflationary or dark energy phase in a AGR cosmology is a semi-global attractor. 
The possibility that scalar tensor gravity could converge to GR has been proposed 
within the context of extended inflation first using numerical techniques [40] and than 
with a more formal proof [41]. The dynamical system approach allows one to see this 
phenomenon in a very clear way, even in the more general case of a non-zero potential. 
It turns out that the nature of the potential, determined in our case by the value of 
n, together with the value of the coupling plays a critical role in the realization of this 
mechanism. 

It is useful to define the deceleration parameter q in terms of the dynamical 
variables: 

9= — -1 = 2* + ^- (33) 

This equation represents a parabola in the phase space that divides the accelerating 
(q < 0) expansion phases from the decelerating (q > 0) ones. Points A and V lie on 
the curve (|33|) as expected by the form of their scale factor solutions. On the other 
hand B and C always lie on the decelerated expansion side of the curve. Instead, for 

£ we have decelerated expansion for 4 - ^i^- <n<2or4<?i<4 + J 'i 2 ^- and 
a decelerated one for < n < 4 - J 12 ^ 1 or n > 4 + W / l2 |~ 1 . 
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Point 



Table 2. Values of the parameter m and the corresponding scalar field solutions 
for the fixed points in vacuum. 

Solutions 



A 
B 



V 



6(l-8C)-8^/3e(12£-l) 
3 



6(l-80+8V3S(12€-l) 



n+2 
2(ra-2) 



h, A 



= Mt-to) 1/2 - m , <^o, A = 
= 0?(i-i o ) 1/2 - m , 



= ± 



4A 
1-12? 



-2 

-2, 2 

4(1- 



12C) 



rt+2 
2(n-2) 



</> = <fi(t-t ) 1 /*, n =-2 

= 0f(i-io) 1 / 2 " m , n^-2, 2, 4 

6Ag 



±2 



^2(1 



12£-1 

2g-l)[3+Q-10)n+2)C] 
Af(n-2)(n-4) 



Table 3. The eigenvalues associated with the fixed points in the vacuum model. 



Point Eigenvalues 



r 2,2] 

4(1 - 120 - V3£(12£-l), 6 - 2(n + 2) + y^WM - 1) 
4(1 - 120 + 8^3C(12C-1), 6 - 2(n + 2) ( 6£ - ^(12^-1) 



(4-n)4-- v /g{3(8-n)ng-4> (4-n)g+ % /g{3(8-n)ng-4> 
(n-2)C ' (n-2)C 

3+(n-10)(ra+2)g 2+2{4+(n-8)n)}£ 
2(n+2)£-l ' 2(n+2)£-l 



^.,2. Asymptotic analysis 

Since the dynamical system (|18j) is not compact, it might admit an asymptotic 
structure that is relevant for the global dynamics. In order to analyze the asymptotic 
features of the phase space, we use the Poincare projection [21,42]. This method 
consists of transforming to the polar coordinates 

i' = fcos'0, y — fsimp (34) 

and setting f = y^- In this way, the asymptotic regime is achieved for r — > 1. Using 
the Poincare projection, the asymptotic form of the dynamical equations Q18p read 



Cosmological dynamics of Scalar-Tensor Gravity 



10 



Table 4. Stability of the fixed points in the vacuum case. The parameters arc 
N ± = 4± ^3(125 - 1)7!, P± = 4±2 V '(125 - l)/3£ and Q± = 4± ^{12$ - 
We use the term 'attractor' to denote a sink, 'repeller' to denote a source and 
'spiral' to denote an attractive spiral. 



Points 




< n < N- 


N- < n < P- 


P- < n< Q- 


A 
B 
C 
V 

e 


saddle 
repeller 
saddle 
spiral 
attractor 


saddle 
repeller 
saddle 
spiral 
repeller 


saddle 
repeller 
repeller 
spiral 
saddle 


saddle 
repeller 
repeller 
attractor 

saddle 




Q- < n < 2 


2 < n < Q+ 


Q+ < n < P+ 


P+ <n< N+ 


A 
B 
C 
V 
£ 


saddle 
repeller 
repeller 

saddle 
attractor 


saddle 
repeller 
repeller 

saddle 
attractor 


saddle 
repeller 
repeller 
attractor 

saddle 


saddle 
repeller 
repeller 
spiral 
saddle 


n> N + 


A 
B 

C 
V 
£ 


saddle 
saddle 

repeller 
spiral 

repeller 









, cos 3 ?/> sin^ 
* = " 6g(l-r)» • (36) 

Note that the radial equation does not depend on the radial coordinate. This means 
that the fixed points can be obtained by only considering the equation for tp' . Setting 
ip' = we obtain the four fixed points listed in Table [5] 

Let us now derive the behaviour of the scale factor at these points. This can be 
done by integrating the Fricdmann equation Q19[) in the above asymptotic limit. The 
points Aoo and C^, correspond to y — > as x — > ±oo respectively. In this limit the x' 
equation (fT5|) reduces to 

x' = ~x\ (37) 

which admits the solution 

x 2 = . (38) 

T - Too 

In the given limit, the evolution equation (|19[) can be written as 
H' 1 , 



which may be solved together with (|38j) to yield the solution 



= [C 1 ±C 2 (t~t )}\ (40) 
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which represents a universe which reaches a maximum radius and then recollapses. 

The points Boo and Poo, instead, correspond to x — > as y — - > ±00. In this limit the 
second of (fT5[) reduces to 

V' = 2( ^-_ 2 l-" »■ , (41) 

which admits the solution 

v = (42 ) 

y 2(l-6e(n-2))(r-roo)' 1 ' 



y, (43) 



Equation ([19]) now reduces to 

H' _ (l-6£(n-2)) 

~H~ ~~ 12f- 1 
and substituting the solution for y, we obtain 

\T-T oo \ = [C 1 ±C 2 {t-t )] 2 . (44) 

The stability of the asymptotic fixed points is summarized in Table O The points 
_4oo and Coo are saddles for every value of n and £. The points Boo and T>oo are 
non-hyperbolic and they can be shown to represent saddle-nodes. This means that 
they behave like saddles or nodes depending on which direction the orbits approach 
them and that a local separatrix exists to divide the different stability domains. In 
our specific case this separatrix corresponds to the equator of the Poincare sphere (i.e. 
our "infinity"), so that effectively Boo behaves as a saddle if n < 2 + l/6£ or n = 2 
and an attractor if n > 2 + l/6£ {n ^ 2), and T>oo behaves like a saddle if n > 2 + l/6£ 
(including n = 2) and an attractor if n < 2 + 1/6£. For both these point the separatrix 
(i.e. our unitary circle) is always attractive so that orbits can bounce off the saddle 
and then approach the point along the unitary circle (see for example Figure [T]) . A 
summary of the stability of the asymptotic fixed points is summarized in Table [51 

The behaviour of the scalar field can be obtained in the same way by substituting 
the solutions for x and y in the appropriate limit of equation (|20p . For the points 
and Coo pnj) we obtain 

<f> = ci cos fcoV^CC* - M) + c 2 sin (co\/3i{t - f )J , (45) 
where cq, c\ and ci are integration constants. For the points Boo and T>oo we obtain 



<f> = Cl COS Co(t — to) 



6g(n-4) 



2(l-6£(n-4)) 



+ C2 sin(c (t-t ) A / 2(1 6 ig ( ^ 4)) , (46) 



for n > 2 + jjj or n > 4, and 



= cie Co( o) V 2(i-es("-4» + C2e Co(t ~* o) V 2(i-6 £ („-4)) ; (47) 
for 2 + ^ < n < 4. 

Since the phase space is two dimensional we can easily draw phase space diagrams 
for the vacuum case. Here wc will limit ourselves to four examples representing the 
global phase space for four specific values of (n, £) (see Figures [TH]) that includes 
the two cases in which the theory admits s a GR attractor. 

|| By "global phase space" we mean the projection of the 8 > n/2 part of the Poincare sphere on the 
plane that contains its equator. 






y 


















(J 


X 









Figure 2. Global phase space 



for n = 2 with £ = 



-1. 
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Table 5. Coordinates, behaviour of the scale factor and stability of the 
asymptotic fixed points in the vacuum model. 

Point ip Behaviour Stability 
-4oo |r - Tool = [Ci ±C 2 (t- t )} 2 saddle 

B °° 2 k - Tool - [Cl ± C 2 (t t Q )] | attractor {n>2+ | } _ {n= 2} 

Coo 7T |t - Tool - [Ci ± C 2 (t - io)] 2 saddle 

o r N ,o f attractor n < 2 + A 

P - T l r_r °°l = [Ci±C 2 (t-t )) 2 [ saddle n>2 + | 



5. The matter case 

As we have seen in the system (j!4j) , the presence of matter in the dynamical system 
equations implies the introduction of another variable: z. Since fi is defined to be 
positive, the definition of the dynamical variables tells us that only negative values of 
z are compatible with attractive gravity. For this reason we will restrict the following 
analysis only to the physically relevant case z ^ 0. 

As in the vacuum case, we can use the constraint (| 15[) to eliminate one of the 
dynamical variables. The system (fl4|) then reduces to: 

x' = [2(1 - 120a + 4(1 - 12£)x 2 - (2 - ^) x 3 + 6£(n - 4)y 

-(1 - 6£(n - 2))xy ~ ±(1 + 3w - 24£)xz + 6£(1 - 3w)z] , 



12£-1 



2(12£ - 1) + (12£ - l)(n - 6)z + 2 ( 2 - A a; 2 



-2(1 - 6£(n - 2))y - (1 + 3w - 24£)z] , (48) 
z' = [(1 + 3w)(12£ - 1) - 6(12$ - l)z 

+2 (2 - x 2 - 2(1 - 6£(n - 2))y - (1 + 3io - 24£)z • 

Since j/' = and z' = are zero for y = and z = 0, the two planes y = and 
2 = corresponds to two invariant submanifolds. The first plane represents classes 
of theories in which the potential is zero, the second one classes of theories for which 
z = and constitutes part of a vacuum invariant submanifold. The structure of the 
total vacuum invariant submanifold can be derived by writing the energy density in 
terms of the dynamical variables [31,32]: 

\i cx zy>*- 2 H ™- 2 . (49) 

It is clear that when z = and y =/= the energy density is zero. However when y = 
and 2^0 the behaviour of (i depends on the value of n: the energy density is zero 
when n > 2, but it is divergent when n < 2. When both y and z are equal to zero, we 
can conclude that [i = only if n > 2. For rt < 2 (|4"9"|) is divergent and /x can only be 
obtained by directly solving the field equations. 
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5.1. Finite analysis 

Setting x' = 0, y' = and z' — we obtain eight fixed points (see Table [S])- The first 
five (.4, B, C, T> and £) sit in the z = plane and have the same (x, y) coordinates of 
the corresponding vacuum fixed points. The coordinates of T and Q are independent 
of ?i and are finite for all £; Ti is a finite fixed point for n ^ and merges with the 
asymptotic fixed subspace £oo for n = 0. 

Merging occurs between the points fi, C, T>, £, Q and Ti. The first four points 
merge in the same way and for the same values of the parameters given in the vacuum 

case. Point H merges with £ for n = ± ( 7 + 3w + / ^ 2 +66^+73)-6(^+i) \ . and ^ 



2 I ' 1 1 V f 
C when n = ^£±lKgl±gzi) 

All the fixed points but A, T>, T are associated with flat solutions. A is associated 
with negative curvature and the sign of the spatial curvature for point D and T 
depends on the value of n and £. 

The stability of the fixed points may be determined using the Hartman-Grobman 
theorem as in the vacuum case. The point A is a saddle for every value of w. Point B 
keeps the same stability found in the vacuum case for w = and 1/3, but it is always 
a saddle when w = 1 . Instead, the stability of points C and T> is the same as in the 
vacuum case for every value of w. Point £ varies its stability with £,n,w and can 
behave like an attractor, a repeller or a saddle as shown in Table [TUJ On the other 
hand, the fixed points T and Q are saddles for every value of the parameters. Finally, 
H is a saddle or a saddle focus for w = and w = 1/3, while for w = 1 it can also be 
a repeller or an anti-spiral. Its stability is summarized in Table [TT1 

As in the vacuum case the coordinates of the fixed points can be used to find 
exact solutions for the evolution of the scale factor. From © and (fTTj) we have 

■ H 2 

H + — = 0, a r , 

01 m 



1 _ O t _ , / l+12g-6fn \ , / l+3^-24g 

1 ix, g^Xj -r 1 x_ 12 £ J Hi T" I 2-24? 



(50) 



0, (3 — 3xi -\~ xf -\ j i_i9f ^ 2/i H ^i_i2f " z * ' (51) 



1-12^ « l 1 1-12? 



(i-<o) 2 

where (xi,yi, Zi) represents the coordinates of the fixed points. For a and f3 different 
from zero, these equations have solutions of the form (|^|) and (I23[) . 

The points ^4, B, C, T> and £ have the same solutions as in the vacuum case, 
since for these points z = 0. In particular the convergence mechanism of scalar tensor 
gravity towards GR is preserved when matter is present. Such a result was expected 
since in [40] and [41] this phenomenon is described when matter is present. The point 
T is associated to the Milne evolution and a constant scalar field 

a r = \ 4>jr = (j) , (52) 

the value of the constant ao being influenced by the parameters £ and w for non-zero 
spatial curvature 

However, a direct check of the equations reveals that this solution is valid only for 
w = —1/3. Since in our model standard matter For point Q we have 

2(8g - 1 + w) _ 3w 2 + 48£w + 16g - 3 

ag ~ 32£ + 3( W 2 -1)' mg ~ 2(3^+32^-3) ' ( } 
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This solution represents a Friedmann like expansion with the exponent depending on 
both £ and w. For < w < 1/3 we have — < otg < g~g^| , for 1/3 < w < 1 we 

have < an < —■ Direct comparison with the cosmological 

3(1 - 8£) - 4 v /3e(12e- 1) " 2 

equations reveals that this solution is physical only for w > 2/3 and w^l, 

Finally, point 7i with n ^ 2 is linked to a solution that a first glance resembles a 

well known Friedmann GR one: 

2rc » + 2 

3(?i — 2){1 + w) 2{n — 2) 

This solution represents a Friedmann-like expansion for n < and n > — - — — — -, a 

3iu + 1 
6(w + 1) 

contraction for < n < 2 and power-law inflation when 2 < n < . When 

3w + 1 

n = 2 we obtain a solution of the form 

a = »o , 4>~4>q, A = 0, (56) 

which corresponds to a static universe with cosmological constant. Also this solution 
represent an effective "GR state" for the scalar tensor cosmology but in this case it is 
not a stable one. This means that in our model of non- vacuum scalar tensor cosmology 
no stable non vacuum GR like solutions are allowed. 

Using (|49|) and the cosmological equations it is easy to conclude that the points 
T and Q both admit vacuum solutions. This is not true for 7i, for which the energy 
density is 

fx = H r^, (57) 
when n ^ 2, where 

4£n 2 



MO = ZH 



_3(n- 2) 2 (1 + w) 2 

and yu and z-h are the y and z coordinates of TL respectively It is clear from the 
expression above that this point does not represent a physical solution for all values 
of n and £ since fio and therefore /i can, in principle, be negative. In order for us to 
determine the values of n and £ for which TL is physical, we have to solve the inequality 
Ho > 0. We obtain 



W +3w )-y ie <n< i( 7 + 3w ) + \j if -( 58 ) 

When n = 2 the cosmological equations reveal that H is instead associated with a 
vacuum (/i = 0) solution. 

As in the vacuum case we can write the deceleration parameter q in terms of the 
dynamical variables via 

1 = -4 - 1 = ^ + eW - ( 1± m^) W - f 1 ^^) *■ (59) 



— 1 ^ 1-125 J y% V 2 - 24 € 

The surface defined above divides the phase space in two volumes representing the 
region in which the expansion is accelerated and the region in which it is decelerated. 
The fixed points A, B, C, T> and £ behave like their vacuum counterparts, since 
they all satisfy z — 0. Point T lies on this surface consistently with the fact that it 
corresponds to Milne evolution. The point Q always lies in the region in which the 
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Table 6. The coordinates and the sign of the spatial curvature of the non vacuum 
fixed points. 



Point Coordinates (x,y, z) K 

A (0, 0, 0) -1 

B ^(-6^-^(12^-1)), 0, o) 

C f2(-6£+ V3£(12£-l)), 0, 0j 

U 1^77=2' 3?(n-2) 2 ' U ) 



£ 



l + [4+n(n-8)]g 
(n-2)^ 



2(n-4)g (12g-l)[3+(r t -10)(n+2)fl 
l-2(n+2)£' 3(l-2(n+2){) 2 
^ | (l+3w) q (12g-l)(l+3m) N \ 



l.o) 



12^ J 16{ 



r / 4g(l-3™) n (12g-l)[3(™-l) 2 -16g(2-3ttQ] \ 

^ I 8£-l+w ' U ' 3(8£-l+™) 2 J U 

o, / 3(l+w) 3-3w 2 -4£[6(l+w)+n(l-3w)] n 
n { n ' iT^Tf ' U 

3(l+M>)+ 2 ?[n 2 -6(l+u))-n(7+3i«)] 
2^1 



Table 7. The exponent a of the scale factor solutions and the energy density for 
the non vacuum case. 

Point a Matter density 

A I n = 

B \ /i = 

3(l-8£)-4y3f(12£-l) 

C -, H = 

3(l-80+V3e(12«-l) 

^ i ' ^ fl - l+[4+n(n-8)]£ ) ^ ~ U 

C 2(n+2)g-l , 2 4 

£ ^ (n-4)(„-2)r /i = 

•T 7 1) ( a = l-16?+ 6 (2-3iu)iu) /i = 

c 2(8g-l+™) n 

y 32£+3(> 2 -l) ' 

^ | 3(n-2Kl+ M ) ' "-T^ 2 | /i = ^ i _ — 

a = ao, n = 2 \/i = 
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Point m 



Table 8. The parameter m and the corresponding scalar field solutions for 
the non vacuum case. The integration constants have been calculated by direct 
substitution in the cosmological equations. 

Solutions 



A 
B 

C 
V 



2(24£+4^3?(12«-l)-3) 

3 

2(-24C+4^/3e(12£-l)+3j 

n+2 



2(n-2) 



h, A = 

6o(t-io) 1/2+m , A = 



= Mt-to) 1/2 ~ 



A = 



= ^(t-<o) 1/2 , » = -2 
= 0f(<-<o) 1/2 - m , 2 



n+2 
2(n-2) 



= 0g(t-to) 1/2 , n=-2 

= 0f(t-t o ) 1/2 - m , 2 







i(3w + 2) 

3w 2 +48g«' + 16£-3 
2(3iu 2 +32{-3) 



(w = -1/3), A = 



A = 0, £ 



3(tti-l)' 
16(3u>-2 



7 -> (w > 2/3, u; ^ 1 ) 



n+2 
2(n-2) 



f = 0*(i-i o ) 1/2 , n = -2 

jM - + o f 24A(1+™) 2 

- =" ^3(l-to 2 )-16C(l+3-(u) 

2(3-4(n+6) + 12(n-2)£m-3w 2 ) ' 
3A(n-2) 2 (l+Mi) 2 



expansion is decelerating. The behaviour of 7i depends on the barotropic factor w 
and the parameter n: it lies in the accelerated evolution region for w = when n < 
and n > 6, for u> = 1/3 when n < and n > 4 and for w = 1 when n < and n > 3. 
Otherwise it lies in the decelerated evolution region. 

5.2. Asymptotic analysis 

We next study the asymptotic behaviour of the system (|48|) using the Poincare 
projection. The compactification of the phase space can be achieved by transforming 
to spherical coordinates 

x = f bvo. 6 cos ij), y = f sin 8 simp, z = f cos 9, (60) 

where f = and f £ [0, oo), 9 £ [0, n] and tp £ [0, 2n]. In the limit r — ► 1 (f — » oo), 
equations Q48p become 

, cos 2 sin 2 



24£ 

cos sin 3 cos 4 "0 



[(cos 20 - 3) sin 2 9 - 4 cos 2 9] , (61) 

(62) 



2 



6£(1 - r) 
. sin 2 9 cos 3 sin ?/i 

^ - m -rr ■ (63) 
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Point Eigenvalues 



A 
B 



V 
£ 
T 

Q 
H 



-2,2,-l-3tul 



4(1-12£)-V3£(12£-1), 6-2(n + 2) 6£ + v /3£(12£-l) , 



3(1 - 80 - ^(12^-1) - 3w 

4(1 - 12^) + 8v/3£(12£-l), 6 - 2(n + 2) (e£ - V3f(12f-1)) , 
3(1 - 80 + v/3?(12C- 1) - 3w 



(4-n)g— y /g{3(8-n)ng-4}' (4-n)g+ - s /g{3(8-n)ng-4}' 6 _„ 



(»-2)« 



(»-2)« 



' n-2 



3iu 



3+Q-10)(n+2)g 2+2{4+(n-8)n)}j 3+2{-6+(n-7)n)}g _ „ ' 
2(n+2)f-l ' 2(n+2)f-l ' 2(n+2)£-l J " ; 



(4+(2-n)(l + 3w)) 



-2g(l-3-i«)-^2£{18£(l+^) 2 -(l-^)(l+3tu) 2 } 



IS, 



-2£(l-3m) + ^2g{18qi+tu) 2 -(l— tu)(l+3m) 2 } 
4C 

16£-(l-m)(l+3t») 3(l-t») 2 -16£(2-3m) 3j> 2 -l)+2g{6+n+3(2-n)tu} 
8{-l+w ' 2(8f-l+«)) ' si-l+w 



-4+(n+2)(l+3m) 3(-m-l)(n+2) + 12-^S(n,^) 3(w-l)(n+2) + 12+- v /g(n,g,^) 
n ' 4n ' 4n 



S(n,£,w) = 3C _1 (-1 + 120 -1 {-36(™ - 1)(1 + TO ) 2 + 3£ [4(1 + w) 2 (12m - 37) 
+n 2 (7 + (2 - 9w)w) + 4n{l + w) (w(19 + 6w) - 17)] + 4£ 2 [324(1 + wf 
+8n 3 (-l + 3w) + 12n(l + w)(29 + 3w) + n 2 (17 - 21^(10 + 3w))] } 



Table 10. Stability of the fixed point £ for the matter case. The parameters 
are N ± = 4 ± ^(12$ = 1)7?. f± = I (7 + 3w) ± ^ S(9^ 2 +66^+73)-6( 1+m ) ^ 
Q± = 4 ± ^(12£ - 







, l-4f 

n < -i-r 


^ < n < N. 


N- < n < P- 


P- 


< n < Q- 


II' 


= 0,1/3 
= 1 


attractor 
attractor 


repeller 
repeller 


saddle 
saddle 




saddle 
saddle 






Q- < n < Q+ 


Q+ <n < P+ 


P+ <n< N+ 




- < n < P\ 


w 
w 


= 0,1/3 
= 1 


attractor 
attractor 


saddle 
saddle 


saddle 
saddle 




repeller 
saddle 



n > P_ 



w = 0,1/3 repeller 
W = 1 repeller 
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Table 11. Stability of the fixed point H of the matter case. The parameters are 
P± = 1(7 + 3w) ± ^J i(9w2 + 66m + 73 > -6(1+^) and Si , which are the real roots of 
the polynomial S(n,£,w) given in Table M for the given value of w. We use the 
term 'anti-spiral' to denote a pure repulsive spiral and 'saddle focus' to denote an 







unstable point w 


ith two complex eig 


cnvalues. 




w = 





n < Si 


Si < n < S2 


S 2 < n < S3 


n > S3 






saddle 

(n/ 3(1 4 -/°: 


saddle focus 


saddle 
(nj4 0,6,P±) 


saddle focus 


w = 


1/3 


n < Si 


Si < n < S2 


n > S 2 








saddle focus 


saddle 
(n#0,4,P±) 


saddle focus 




W = 


1 


n < Si 


Si < n< P+ 


P+ < n < S 2 


S 2 < n < S 3 






saddle focus 


saddle 
(n/0,3 ) P_) 


repeller 


anti-spiral 






S 3 < n < 6 


n > 6 










repeller 


saddle 







As in the vacuum case the radial equation does not contain the radial coordinate, 
so that the fixed points can be obtained using only the angular equations. Setting 
tp' = and 9' = 0, we obtain four fixed points and a fixed subspace which are listed 
in Table QJ] 

The points and Boo he on the plane z = and therefore their solutions are 
the same as the vacuum points Aoo and Coo , respectively. Points Coo and T>oo represent 
the poles of the Poincare sphere and it is easy to prove that they are linked to the 
same solutions at Aoo and Boo- The solution for the scalar field instead reads 



Co(t-t )\ — 9_?4<r , c (t-t )-' 



= Cie -v- -V 2-24? +C2 e '- 2-24C _ (64) 

In addition to the fixed points above, we found two fixed subspaces L\ and £2, 
that contains all the points with -0 = it/2 and -0 = 3ir/2, respectively. These subspaces 
corresponds to x — > 0, y/z ±tan#n, respectively. The equations for y' and z' in 
this limit reduce to 

y = [2(1 " 6f(n " 2)) ± (1 + 3w ~ 24?) cot 9o] ' (65) 

[(l + 3w-24f)±2(l-6£(n-2))tan0o] , (66) 



1 - 12£ 
which admit the solutions 



= m-i 

y (r -r^) [2(1 -6£(n- 2)) ± (1 + 3w - 24£) cot 0„] , 1 j 



12£- 1 



(r - r^) [(1 + 3w - 24£) ± 2(1 - 6£(n - 2)) tan0 o ] 
In the limit above, (|5D|) becomes 



(68) 



H = - 
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which implies 

\T-T oo \ = [C 1 ±C 2 (t-t )} 2 . (70) 

The solution for the scalar field in the fixed subspaccs £1 and £2, can be found in 
similar way. In the limit above, using (|67|) . (|68|) and (|70jl in (fSTj) we find 

I + A C 2 = , (71) 

where 

= 6g(l - 3uQ cos fl ± 12g(n - 4) sin 9 

2(l + 3w-240cos6> ±4(l-6^(n-2))sin6'o ' 1 j 

We find the following solutions 

ci cos{^/\A\c (t - t )) + c 2 sm(^/\A\c {t - < )), if A < 0, 



ci 



e \/3co(t-to) _j_ ^gx/Acott-to)^ ifA>0. 



(73) 



The stability of these asymptotic fixed points can be deduced by analyzing the 
stability with respect to the angular coordinates, and from the sign of r' . The points 
„4oo and B^, are stable nodes for all values of n and £. Points and 2?oo instead are 
always saddles. 

All the other points including the fixed subspaccs C\ and £2, have both 
eigenvalues equal to zero. In order to derive the stability of these points we have 
to analyze the effect of the non-linear contributions have on dynamical equations. 
This can be done by Taylor developing the R.H.S of the dynamical equations around 
the fixed point up to the first non-zero order and then directly the resulting system. 

For the points and 2?oo we obtain the solutions: 

ih c = c 2 , 9c = — , (74) 

r sin (^0) cos 3 (ip ) + 6£ci ' 

1pD x = C 2 , 9 Daa = 7T . . (75) 

rsin (tpo) cos 3 (yj ) + 6£ c\ 

This result tells us that when we choose initial conditions around these points, the 
evolution of the universe will follow an orbit with constant ip, with r increasing 
and 9 approaching or n. If we also consider the behaviour of the radial equation we 
conclude that these points behave like saddles. 
For the fixed subspaces £1 and £2, we find 



V r cos 6» sin 3 9 + 96^ci 2 

9(t) =0 O - ^|^ + C2[t (sin 40o-2 sin20 o )- 768eci)]* csc2e ° . (77) 

It is clear that, for r increasing, the solution for approaches ir/2 for £1 and 37r/2 for 
£2, while 9 increases. The radial behaviour is very complicated and depends critically 
on the value of the coordinate 9 , w, n and £ (we will not show this dependence here). 
We may however conclude that this subspace is a saddle or an attractor depending on 
the value of 9 . 



% This behaviour is reversed in case of a contracting cosmology because in this case r effectively 
changes sign when a is decreasing. 
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Table 12. Coordinates, behaviours and stability of the asymptotic fixed points 
in the non vacuum case. The value of the parameter A is given in H72I I. 



Point 


(M) 


Scale factor 




Stability 


Boo 
Coo 


(|,o) 
(M 

(0,0) 

(?T, 7r) 


|r-Too| = [Ci±C a (t- 
|t-Too| = [Ci±0,(t- 
Ir-Tool = [Ci±C a (t- 
|r-Too| = [Ci±C 2 (t- 


*o)] 2 
*o)] 2 
io)] 2 
*o)f 


attractor 
attractor 
saddle 
saddle 


A 

£2 


(00, f ) 
(^O, 3 1 ) 


k-Tool = [d±C 2 (i- 
|r-Too| = [Ci±(7 2 (t- 


io)] 2 
*o)] 2 





The phase space in the matter case is 3-dimensional and therefore cannot be 
visualized as easily as its vacuum counterpart. For this reasons we will not give here 
any sketch of the phase space and we refer the reader to the next section for an analysis 
of the results derived above. 

6. Discussion and conclusions 

In this paper, the dynamics of STG FLRW cosmological models has been studied using 
a phase space analysis. We have considered a generic non-minimally coupled theory 
of gravity where the coupling and the potential are powers of the scalar field, both in 
a vacuum and in presence of a perfect fluid. The set of parameters characterizing the 
cosmological models are {£, A, n, to}, i.e. the coupling constant, a constant parameter, 
the power of the self-interacting potential and the barotropic index of the perfect 
matter fluid, respectively. The phase-space is 3-dimensional in absence of matter while 
it is 4-dimensional in presence of matter, but in both cases the FLRW Hamiltonian 
constraint allows one to reduce their dimensionality. Our investigation considered 
the existence and local stability of critical points (finite analysis) and the asymptotic 
analysis via the Poincare projection. 

We identified five finite fixed points in the vacuum case and eight in the matter 
case. In the vacuum case, there are four asymptotic stability points corresponding 
to the four intersections of the axes with the unitary Poincare circle. In the matter 
case, we have to consider a unitary 3-sphere and in additin to the vacuum asymptotic 
points we find two more points and two "fixed subspaces" . The stability of fixed points 
strictly depends on the values of the above parameters and in particular on £ and n, 
i.e. the coupling and the power of self-interacting potential. 

In the vacuum case, the two dimensional phase space is divided in two halves by 
the invariant submanifold associated with V(</>) = 0. Of the nine fixed points we found, 
five of them (A, B,C,Aoo, Coo) ar e permanently on this invariant submanifold, two of 
them Boo, Coo do not change their position and the other two "move" in the phase 
plane depending on the values of the parameters. In particular, T> is characterized by 
q = and its position is on the curve described by the equation given in (|33[) and 
£ is associated with a flat spatial geometry and always lies on the curve described 
by the constraint ([17]) with K = 0. Even if matter is not present the scalar field is 
able to induce an expanding Friedmann-like evolution. However, the value of a in 
(|26p and (|2~T|) reveals that at these fixed points the scale factor cannot grow faster 
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than i 1 / 2 . In addition, these solutions are unstable for every value of the parameters. 
The point £ admits the widest spectrum of behaviours. Depending on the values of 
n and £ its solution can represent an inflationary phase, a Friedmann-like phase or a 
contraction. A comparison with the stability analysis reveals that £ is an attractor 
only when it corresponds to contracting solutions or power law inflation. This has 
two consequences: (i) in the scalar field dominated regime our model of scalar tensor 
cosmology admits an inflationary phase as an attractor even if this attractor is not 
global and not unique; (ii) since there is no value of the parameters for which the £ 
is a saddle, this model does not admit a transient inflationary phase and cannot solve 
the graceful exit problem. 

The case n = 4 is particularly interesting because, for this value of n, £ is 
associated with a GR (cf> = const.) de Sitter solution and is an attractor. This fact 
has two main consequences: (i) it shows that with a non-minimally coupled scalar 
field with quartic potential, the early time cosmology evolves towards an inflationary 
phase which is indistinguishable from a GR one; and (ii) it gives us an independent 
confirmation of the idea that a scalar tensor theory of gravity can evolve towards GR. 
The difference is that when a potential is present, the realization of such mechanism 
is strictly related to the form of the potential. 

In the matter case, we found three new finite fixed points which are physical 
only for specific values of the parameters. In particular, T is never physical because 
it satisfies the cosmological equations only for a negative w. Point Q represents a 
physical vacuum solution only for w > 2/3 and w =/= I. For these values of w, 
this point is associated with an expansion whose rate cannot be higher than the 
radiation dominated GR-FLRW solution. Point H is associated to a non-vacuum 
solution which resembles a well known Fricdmann-GR solution. As n varies, this 
solution can represent power law inflation, a decelerated Friedmann solution and a 
contraction. For n = 2, this point is linked to a GR-like state, but since it is unstable 
for every value of the parameters, we conclude that there is no way for this class 
of scalar tensor cosmologies to approach a stable non-vacuum GR state. In fact, a 
quick look to the table above reveals that there is no value of the parameters for 
which any of the finite fixed points are stable. This means that the only non-vacuum 
attractors for this class of theories are the asymptotic fixed points and their associated 
Lemaitre solution. Therefore we can conclude that the class of cosmologies we treated 
are doomed either to approach an effectively vacuum state (probably corresponding 
to thermal death) that corresponds to one of the vacuum attractors or to recollapses 
towards a Big Crunch. 

The most interesting orbits are the ones that "travel" between point £ and 
TL. This is because they could represent cosmic histories in which a Friedmann- like 
cosmology enters naturally in a phase of accelerated expansion or cosmic histories 
in which an unstable inflationary phase is followed by a second inflationary phase. 
The first type of orbits is interesting because they can in principle help solve the 
incompatibility between the evolution towards a Dark Energy era and the formation 
of large scale structure. The second ones are interesting because they potentially unify 
"dark" scalar fields and the inflaton within a single scheme. 

Since the phase space is three dimensional it is not easy to check if such orbits 
actually exist without the use of numerical techniques. However, our results allow one 
to give some necessary condition for these orbits to exist and to rule out some of them. 
For example, it is easy to see that there is no value of the parameters for which £ can 
represent an unstable Friedmann solution and TL an inflationary phase. On the other 
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hand, for 4 < n < Q + and w = 0, 1/3 we have that TL corresponds to an unstable 
decelerated expansion, S correspond to a stable power law inflation and these points 
are not separated by any invariant submanifold. This means that, in principle cosmic 
histories exist for which a transient Friedmann evolution approaches to a power law 
inflationary phase in a natural way. 

If 4 < n < 6(1 + w)/(3w + 1) and w = 0, 1/3 another interesting set of cosmic 
histories is possible in which we have two inflationary phases, a first one which is 
unstable associated with 7i and a second one which is stable associated with £. 
Although only a detailed analysis of these cosmic histories can reveal if this last group 
of cosmic histories also include the deceleration phase necessary for the realization 
of standard cosmology, this scenario is interesting because it shows that in STG the 
non-minimally coupled scalar field can act as both the inflaton and dark energy. Such 
a behavior has been also found in a class of higher order gravity models [43] and 
definitely deserves a more careful investigation. Such a study, will be the topic of a a 
forthcoming paper. 

As final comment, it is worth stressing the connection between the model 
presented above and the String-Dilaton action. Using the transformation ([5]), which 
is only a reparametrization for the scalar field, we can pass from the action (p} 
to ((6|). This means that all the exact solutions for the scale factor that we have 
derived together with their stability are also solutions for the String-Dilaton action. 
Therefore, our analysis also allows us to give details of the dynamics of low-energy 
string cosmology which could be of great interest in the quest for finding observational 
constraints for this theory. 

Acknowledgments 

This research was supported by the National Research Foundation (South Africa) 
and the Italian Ministero Degli Affari Esteri-DG per la Promozione e Cooperazione 
Culturale under the joint Italy/ South Africa Science and Technology agreement. The 
authors wish to thank A Starobinsky and G Vcnturi for their useful comments. 

References 

[1] Weinberg S 1972 Gravitation and Cosmology (New York: Wiley) 

[2] Fujii Y and Maeda K 2003 The Scalar-Tensor Theory of Gravity (Cambridge: Cambridge 
University Press) 

[3] Jordan P 1955 Schwerkraft und Weltall (Braunschwig: Frierich Vieweg und sohn) 
[4] Brans C and Dicke R H 1961 Phys. Rev. 124 925 

[5] S. Capozziello, S. Nojiri, S. D. Odintsov and A. Troisi, Phys. Lett. B 639 (2006) 135 
[arXiv:astro-ph/060443"T]; S. Cap ozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 634 
(2006) 93 |arXiv:hep-th/0512118] . 

[6] Buchbinder I L, Odintsov S D, and Shapiro I L 1992 Effective Action in Quantum Gravity 
(Bristol: IOP Publishing) 

[7] Guth A 1981 Phys. Rev. D 23 347 

[8] La D and Steinhardt P J 1989 Phys. Rev. Lett. 62 376 
[9] La D, Steinhardt P J and Bertschinger E 1989 Phys. Lett. B231 231 
[10] Peebles P J E and Rathra B 2003 Rev. Mod. Phys. 75 559 

[11] Cardone V F, Troisi A, and Capozziello S 2005 Phys. Rev. D72 043501 and references therein 

[12] Capozziello S and do Ritis R 1993 Int. Jou. Mod. Phys. D 2 367 

[13] Capozziello S et al 1996 La Riv. del N. Cim. 4 1 

[14] Capozziello S and de Ritis R 1994 Class. Quantum Grav. 11 107 

[15] For a review on the dynamical system theory a useful book is G Sansone and R Conti 1956 
"Equazioni Diffcrcnziali Non Lincari" (Edizioni Cremonese, Roma) ; for the application of 



Cosmological dynamics of Scalar-Tensor Gravity 



25 



dynamical system to cosmology see Wainwright J and Ellis G F R (ed) 1997 Dynamical 
systems in cosmology (Cambridge: Cambridge University Press) and references therein 

[16] Qualitative analysis of cosmological models in Brans-Dicke theory, C.Romero, H.P. Oliveira and 
J.R.T. Mello Neto, Astrophysics and Space Science 158, 229 (1989) 

[17] Kolitch S J and Eardley D M 1995 Ann. Phys. (NY) 241 128 

[18] Kolitch S J 1996 Ann. Phys. (NY) 246 121 

[19] Santos C and Gregory R 1997 Ann. Phys. (NY) 258 111 

[20] Holden D J and Wands D 1998 Class. Quantum Grav. 15 3271 

[21] Amendola L, Litterio M and Occhionero F 1990 Int. Journ. Mod. Phys. A 5 3861 

[22] D. J. Holden and D. Wands, Phys. Rev. D 61 (2000) 043506 |arXiv :gr-qc/9908026 . 

[23] L. Amendola, Phys. Rev. D 60 (1999) 043501 [arXiv:astro-ph/9904120f 

[24] J. D. Barrow and K. i. Maeda, Nucl. Phys. B 341 (1990) 294. 

[25] R. Gannouji, D. Polarski , A. Ranquet and A. A. Starobinsky, JCAP 0609 (2006) 016 

[arXiv: astro-ph /0606287] . 

[26] Foster S 1998 Class. Quantum Grav. 15 3485 

[27] Billyard A, Colcy A A, Ibanez J 1999 Phys. Rev. D59 023507 

[28] Gunzig E et al. 2000 Class. Quantum Grav. 17 1783 

[29] Gunzig E et al. 2001 Phys. Rev. D63 067301 

[30] Saa A et al. 2001 Int. Journ. Theor. Phys. 40 2295 

[31] Carloni S, Dunsby P K S, Capozziello S and Troisi A 2005 Class. Quantum Grav. 22 4839 

[32] Leach J A, Carloni S and Dunsby P K S 2006 Class. Quantum Grav. 23 4915 

[33] Jarv L, Kuusk P, M Saal 2007 Phys. Rev. D75 023505 

[34] Faraoni V 2005 Ann. Phys. (NY) 317 366 

[35] Capozziello S, de Ritis R, and Rubano C 1993 Phys. Lett. A 177 8 

[36] Birrell N D and Davies P C W 1982 Quantum Fields in Curved Space (Cambridge: Cambridge 
University Press) 

[37] Capozziello S, de Ritis R, Marino A A 1997 Class. Quantum Grav. 14 3243 

[38] Leach J A, PhD thesis 2008 

[39] Hartman P 1964 Ordinary Differential Equations (New York: Wiley) 

[40] J. Garcia-Bellido and M. Quiros, Phys. Lett. B 243 (1990) 45. 

[41] T. Damour and K. Nordtvedt, Phys. Rev. D 48 (1993) 3436. 

[42] Capozziello S, Occhionero F and Amendola L 1993 Int. Journ. Mod. Phys. D 1 615 

[43] M. Abdelwahab, S. Carloni and P. K. S. Dunsby, larXiv:0706.1375l fer-qcl . 



